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Abstract 

We consider the ground state properties of an inhomogeneous two-dimensional Bose gas with 
a repulsive, short range pair interaction and an external confining potential. In the limit when 
the particle number N is large but pa 2 is small, where p is the average particle density and a the 
scattering length, the ground state energy and density are rigorously shown to be given to leading 
order by a Gross-Pitaevskii (GP) energy functional with a coupling constant g ~ 1/| ln(pa 2 )|. 
In contrast to the 3D case the coupling constant depends on N through the mean density. The 
GP energy per particle depends only on Ng. In 2D this parameter is typically so large that 
the gradient term in the GP energy functional is negligible and the simpler description by a 
Thomas-Fermi type functional is adequate. 

1 Introduction 

Motivated by recent experimental realizations of Bose-Einstein condensation the theory of dilute, 
inhomogeneous Bose gases is currently a subject of intensive studies. Most of this work is based on 
the assumption that the ground state properties are well described by the Gross-Pitaevskii (GP) 
energy functional (see the review article Q ) . A rigorous derivation of this functional from the basic 
many-body Hamiltonian in an appropriate limit is not a simple matter, however, and has only been 
achieved recently for bosons with a short range, repulsive interaction in three spatial dimensions 

The present paper is concerned with the justification of the GP functional in two spatial di- 
mensions. Several new issues arise. One is the form of the nonlinear interaction term in the energy 
functional for the GP wave function <J>. In three dimensions this term is Ana J |$| 4 , where a is 
the scattering length of the interaction potential. The rationale is the well known formula for 
the energy density of a homogeneous Bose gas, which, for dilute gases with particle density p, is 
Anap 2 . This fact has been 'known' since the early 50's but a rigorous proof is fairly recent Q. In 
two dimensions the corresponding formula is 47rp 2 | \n{pa 2 )\~ l as proved in [Q] by extension of the 
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method of p]. The formula was first stated by Schick [||; other early references to this formula 
arc |, @, H- It would seem natural to consider Air J |<3?| 4 | ln(|<3?| 2 a 2 )| 1 as the interaction 



term in the GP functional, and this has indeed been suggested in |TT], [PJ] . Such a term, however, is 
unnecessarily complicated for the purpose of leading order calculations. In fact, since the logarithm 
varies only slowly it turns out that one can use the same form as in the three dimensional case, 
but with an appropriate dimensionless coupling constant g replacing the scattering length, and still 
retain an exact theory (to leading order in p). 

It is often assumed that a justification of the GP functional depends on the existence of Bose 
Einstein condensation. Several remarks can be made about this: 1. We neither assume nor prove 
the existence of BE condensation, but we do demonstrate a kind of condensation over a distance 
that is fixed (i.e., non-thermodynamic) but whose length goes to infinity as the density goes to 
zero; 2. BE condensation does not exist in two dimensions when the temperature is positive, but it 
can, and most likely does, exist in the ground state; 3. In any event, when the density is low and 
the temperature is zero it appears to be likely that the system can be described for many purposes 
in terms of only a few macroscopic order parameters such as the density and phase - at least this 
is true for the dependence of the ground state energy and density upon an external potential. 

The functional we shall consider is 

£ GP {$} = J (|V$(x)| 2 + V(x)|$(x)| 2 + 47ra|$(x)| 4 )d 2 x, (1.1) 

where V is the external confining potential and all integrals are over M 2 . 

The choice of g is an issue on which there has not been unanimous opinion in the recent papers 
l|, ||, 0, P, ||, ||, HH on this subject. We shall prove that a right choice is g = | In (pa 2 ) | 1 where 



p is a mean density that will be defined more precisely below. This mean density depends on the 
particle number N, which implies that the scaling properties of the GP functional are quite different 
in two and three dimensions. In the three-dimensional case the natural parameter is Na/a osc , with 
a osc being the length scale defined by the external confining potential. If a/a osc is scaled like 1/N 
as N — > oo this parameter is fixed and the gradient term J \ S/§\ 2 in the GP functional is of the 
same order as the other terms. In two dimensions the corresponding parameter is N\ ln(pa 2 )\ . 
For a quadratic external potential p behaves like A^ 1 / 2 /a 2 sc and hence the parameter can only be 
kept fixed if a/a osc decreases exponentially with N. A slower decrease means that the parameter 
tends to infinity. This corresponds to the so-called Thomas Fermi (TF) limit where the gradient 
term has been dropped altogether and the functional is 

£ TF [P\ = j (V(x)p(x) + 4vrop(x) 2 ) d 2 x, (1.2) 



defined for nonnegative functions p. Our main result, stated in Theorems |l.3| and [Lj below, is 



that minimization of (1.2) reproduces correctly the ground state energy and density of the many- 
body Hamiltonian in the limit when ./V — ► oo, pa 2 — > 0, but N\ ln(/?a 2 )| _1 — > oo. Only in the 
exceptional situation that N\ ln(pa 2 )| _1 stays bounded is there need for the full GP functional 
([El]), cf. Theorems [EJ and |CJ, 



We shall now describe the setting more precisely. The starting point is the Hamiltonian for N 
identical bosons in an external potential V and with pair interaction v, 

N 

H(N) = E (" v * + + 2>( x * - x ;)< (i-^) 

i=l i<j 

acting on the totally symmetric wave functions in <g> N L 2 (M 2 ). Units have been chosen so that 
h = 2m = 1, where m is the particle mass. We assume that v is nonnegative and spherically 



LSY 19/10/00 



3 



symmetric with a finite scattering length a. (For the definition of scattering length in two dimensions 
see the appendix.) The external potential should be continuous and tend to oo as |x| — > oo. It is 
then possible and convenient to shift the energy scale so that min x V(x) = 0. For the TF limit 
theorem we shall require some additional properties of V to be specified later. 

The ground state energy u> of the one-particle operator — V 2 + V is a natural energy unit and 
gives rise to the length unit a osc = u ' . In the sequel we shall be considering a limit where cl j ct osc 
tends to zero while TV — > oo. Experimentally a/a osc can be changed in two ways: One can either 
vary a osc or a. The first alternative is usually simpler in practice but very recently a direct tuning of 
the scattering length itself has also been shown to be feasible Mathematically, both alternatives 
are equivalent, of course. The first corresponds to writing V(x) = a~ 2 c V(x/a osc ) and keeping V and 
v fixed. The second corresponds to writing the interaction potential as v (x) = a~ 2 t)(x/a), where v 
has unit scattering length, and keeping V and v fixed. This is equivalent to the first, since for given 
V and v the ground state energy of (|Q|), measured in units of u, depends only on N and a/a osc . 
In the dilute limit when o is much smaller than the mean particle distance, the energy becomes 
independent of v. 

We shall measure all energies in terms of to and lengths in terms of a osc and regard V and v as 
fixed. The notation E^ M (N,a) for the ground state energy of (|l.3j ) is then justified. 
The quantum mechanical particle density is defined by 

pf l a ( x )=Nj |^)(x,x 2 ,...,x JV )| 2 d 2 x 2 ...d 2 x iV , (1.4) 

where vpW is a ground state for ( |l,3j ). 

The GP functional fll-l| ) has an obvious domain of definition (cf. Eq. (2.1) in ||). The infimum 
of f GP [$] under the condition f |<1>| 2 = N will be denoted by E GP (N,g). The infimum is obtained 
for a unique, positive function, denoted ^% P g , and the GP density is defined as p GP g (x) = (x) 2 . 

The ground state energy of the TF functional ( |1.2|) with the subsidiary condition J p = N is 
denoted E TF (N,g). The corresponding minimizer can be written explicitly; it is 

P^(x) = ^-[/. TF -y(x)] + , (i.5) 

where [t]+ = max{i, 0} and p TF is chosen so that the normalization condition J pjf g = N holds. 

We now define the mean density p as the average of the TF density pj^ l at coupling constant 
g = 1, weighted with N~ 1 pjf l , i.e., 



P 



1 j p^(x) 2 d 2 x. (1.6) 



It is clear that p depends on N and when we wish to emphasize this we write pN- The definition 
( [LCD has the advantage that p is easily computed; for instance, if V(x) ~ |x| s for some s > 0, then 
Pn ~ N s /( s+2 \ It may appear more natural to define p self-consistently as p = jj J pj^ g (x) 2 d 2 x 
with g = | ln(pa 2 )| _1 , which amounts to solving a nonlinear equation for p. Also, the TF density 
could be replaced by the GP density. However, since p will only appear under a logarithm such 
sophisticated definitions are not needed for the leading order result we are after. The simple formula 
(Oj) is adequate for our purpose, but it should be kept in mind that the self-consistent definition 
may be relevant in computations beyond the leading order. 

With this notation we can now state the two dimensional analogue of Theorem 1.1 in Q. 

Theorem 1.1 (GP limit for the energy). //, for N — ► oo, o?pN — > with N/\ ln(a 2 /?7v)| fixed, 
then 

E QM (N,a) 

N^L EGP(N,l/\ln(a 2 p N )\) = L (L?) 
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The corresponding theorem for the density, c.f. Theorem 1.2 in H, is 

Theorem 1.2 (GP limit for the density). //, for N — > oo, a 2 p~N — > with 7 = N/\ ln(a 2 pjv)| 
fixed, then 

Jjm ^P$J(x) = ^(x) (1.8) 

in the sense of weak convergence in L 1 (IR 2 ). 

These theorems, however, are not particularly useful in the two dimensional case, because the 
hypothesis that N/\ln(a 2 pm)\ stays bounded requires an exponential decrease of a with N. As 
remarked above, the TF limit, where N/\ ln(a 2 /9/v)| ~~ * °°; is much more relevant. Our treatment of 
this limit requires that V is asymptotically homogeneous and sufficiently regular in a sense made 
precise below. This condition can be relaxed, but it seems adequate for most practical applications 
and simplifies things considerably. 

Definition 1.1. We say that V is asymptotically homogeneous of order s > if there is a function 
W with W(x) / for x ^ such that 

A-*y(Ax) - W(x) 

l + \W(x)\ >° aS A ^°° (L9) 

and the convergence is uniform in x. 

The function W is clearly uniquely determined and homogeneous of order s, i.e., W(Xx) = 
X s W(x) for all A > 0. 

Theorem 1.3 (TF limit for the energy). Suppose V is asymptotically homogeneous of order 
s > and its scaling limit W is locally Holder continuous, i.e., |W(x) — W(y)| < (const.) |x — y\ a 
for |x|, |y| = 1 for some fixed a > 0. //, for N — > oo, a 2 p~N — ► but N/\ ln(a 2 p~N)\ 00, then 

E QM (N,a) 

lim . 7— ' n = 1. (1-10) 

N^ooE TF (N,l/\ln(a 2 p N )\) V ; 

To state the corresponding theorem for the density we need the minimizer of ( |1.2| ) with g = 1, 
V replaced by W, and normalization J p = 1. We shall denote this minimizer bypf\; an explicit 
formula is 

^J(x) = i-[/i TF -Ty(x)] + , (i.ii) 

where /i TF is determined by the normalization condition. 



Theorem 1.4 (TF limit for the density). LetV satisfy the same hypothesis as in Theorem 1.5 . 
If, for N — > oo, a 2 pjv — ► 6ui 7 = N/\ ln(a 2 p^)\ — > 00, t/ien 

lim ^^pK(7 1/{s+2) x)=^J(x) (1.12) 
in i/ie sense of weak convergence in L 1 (M 2 ). 

Remark 1.1. For large iV, pjy behaves like (const. )jV S// ( s+2 ). Moreover, prefactors are unimportant 



in the limit N — > 00, because p^v stands under a logarithm. Hence Theorems |L3| and |1.4j could also 
be stated with N s /( s+2 *> in place of pN- 

The proofs of these theorems follow from upper and lower bounds on the ground state energy 
E^ M (N, a) that are derived in Sections || and ||. For these bounds some properties of the minimizers 
of the functionals (|1 . 1|) and (|l,2j), discussed in the following section, are needed. 
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2 GP and TF theory 

In this section we consider the functionals (|Ll]) and ( |1.2| ) with an arbitrary positive coupling 
constant g. Existence and uniqueness of minimizers is shown in the same way as in Theorem II. 1 
in m. The GP energy E GP (N,g) has the simple scaling property E GP (N,g) = NE GP (l,Ng). 
Likewise, iV _1//2 <3?^ = <p GP depends only on 

7 = Ng (2.1) 

and satisfies the normalization condition f |<^> GP | 2 = 1. The variational equation (GP equation) for 
the GP minimization problem, written in terms of <f> GP , is 



-A0f + + 8vr 7 (0f) 3 = m OP (7)^ P , (2-2) 

.GP/ 



where the Lagrange multiplier (chemical potential) fi (7) is determined by the subsidiary nor- 
malization condition. Multiplying (|2.2| ) with GP and integrating we obtain 

M GP (7) = E GP (l, 7) + 4^7 / ^ GP (x) 4 d 2 x. (2.3) 

For the upper bound on the quantum mechanical energy in the next section we shall need a 
bound on the absolute value of the minimizer 4> GP . 



Lemma 2.1 (Upper bound for the GP minimizer). 

Proof. 4> GP is a continuous and positive function that satisfies the variational equation 

-A^P + u^r = ^gp (2 5) 

with U = V + 8vr 7 (0 GP ) 2 . Le t B = {x|0 GP (x) 2 > /j GP /(87r 7 )}. Since V > we see that 
-A<p GP < on B, i.e., 4> GP is subharmonic on B. Hence <j) GP achieves its maximum on the 
boundary of B, where GP (x) 2 = // GP /(87r 7 ), so B is empty. □ 



The ground state energy E TF (N,g) of the TF functional (1.2) scales in the same way as 
E GP (N,g), i.e., E' TF (N,g) = NE TF (1, Ng), and the corresponding minimizer pjf is equal to 
NpJJjg. For short, we shall denote pf F by pJj F . By ([Dj) we have 



P 7(x) = ^-[/, TF ( 7 )-y(x)] + , (2.6) 
with the chemical potential // TF ( 7 ) determined by the normalization condition J pI F = 1. In the 



same way as in (|2.3|) we have 

^ F ( 7 ) = £ TF (1, 7 ) + 4^ 7 J p FF (x) 2 d 2 x. (2.7) 

The chemical potential can also be computed from a variational principle: 
Lemma 2.2 (Variational principle for p TF ). 

^ TF ( 7 )= inf /Vp + 87r 7 ||H|oo (2.8) 
p>o,J p=i J 
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Proof. Obviously, the infimum is achieved for a multiple of a characteristic function for some 
measurable set 1Z C R 2 . If \1Z\ denotes the Lebesgue measure of 1Z, then 



inf / V p + 87r7||p| 
fp=iJ 



= i i i (L v+8 ^)w\ (2 - 9) 

= inf (J (V - ^ TF ( 7 )) + 8vr 7 + /, TF ( 7 )!^|) (2.10) 

Now f n {V — Ai TF (7)) > — 87T7, with equality for 

(x|y(x) < ^ TF ( 7 )} C U C (x|y(x) < ^ TF ( 7 )} • (2.11) 

□ 

Corollary 2.1 (Properties of // TF (7)). /i TF (7) is a concave and monotonously increasing func- 
tion of "<( with p TF (0) = 0. Hence A iTF (7)/7 is decreasing in 7. Moreover, /i TF (7) — > 00 and 
M TF (t)/7 ^ is 7 ^ 00. 



Proof. Immediate consequences of Lemma 2.2 , using that min x l/(x) = and limixi^oo V(x) = 
00. □ 

Note that since E TF (l,j) > ^ TF (7) we also see that .E TF (1,7) — > 00 with 7. In this limit the 
GP energy converges to the TF energy, provided the external potential satisfies a mild regularity 
and growth condition: 

Lemma 2.3 (TF limit of the GP energy). Suppose for some constants a > 0, L\ and L2 

\V(x) - V(y)\ < Li|x - y\ a e L ^-^(l + V(x)). (2.12) 

Then 

E GP (1 7) 

lim ) ' r = 1. (2.13) 

Proof. It is clear that E TF (l,^f) < E GP (1,7). For the other direction, we use (j e */a rF ) 1//2 as a test 
function for £ GP , where 

jt(x)= 2i exp H |x| )- (2 - i4) 

Note that J j e = 1 and |Vj e | = e je- Therefore 

where we have used convexity for the last term. Moreover, 

J (Je * V - V)p™ = JI d 2 xd 2 yj £ (x - y) (V(x) - V(y)) p™(x) 

-2^2/ / d2xd2 y|x-y| a e(-^ 1+L2 )i x -yi(i + y(x))^F(x) 

< (const.) e° (l + £ TF (l, 7 )) , (2.16) 
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as long as e < L 2 1 . So we have 

£ GP (1,7) < (1 + (const.) e")^ (1, 7 ) + ^ + (const.) e°. (2.17) 
Optimizing over e gives as a final result 

£ GP (1, 7 ) < £ TF (1, 7 ) (l + (const.)^ TF (l, 7 )- Q /( a+2 )) . (2.18) 

□ 

Condition ( [2.12| ) is in particular fulfilled if V is homogeneous of some order s > and locally 
Holder continuous. In this case, 

£TF (1)7)=7 s/(s+2)£TF (M) (219) 



and 



By fl2.7|) we also have 



7 2/( S+2 ) p TF (7 l/(s + 2) x) = p TF (x) _ (2 _ 20) 



M TF (7)=7 V( S+ 2) M TF (1) _ (2 _ 21) 



If V is asymptotically homogeneous with a locally Holder continuous limiting function W, we 
can prove corresponding formulas for the limit 7 — ► 00. This is the content of the next theorem, 
where we have included results on the GP — > TF limit as well: 

Theorem 2.1 (Scaling limits). Suppose V satisfies the condition of Theorem 1.5. Let .E TF (1, 1) 
be the minimum of the TF functional ( [j.^ j with g = 1 and N = 1 and 1/ replaced by W, and let 
Pi\ be the corresponding minimizer. Then 

(i) lim 7 ^ 00 £; GP (l 5 7)/7 s /(^2) =lim 7 ^ OOJ E;TF (1) = £TF (M) _ 

(ii) lim 7 ^ 0O 7 2 /( s+2 )p G P(7 1 /(«+2) x ) = ^J(x), strongly in L 2 (R 2 ). 
(hi) lim 7 ^oo 7 2 /( s+2 ) / a rF (7 1 /( s+2 )x) = pj\(x), uniformly in x. 

Proof. With the demanded properties of V, ( p.!3| ) holds. Using this and (1.9) one easily verifies 
(i). Moreover, 7 2// ^ +2 )p GF (7 1 / ( ' s+2 ^x) is a minimizing sequence for the functional in question, so 
we can conclude as in Theorem II. 2 in [Q] that it converges to p^ F (x) strongly in L 2 , proving (ii). 
(Remark: In Eq. (2.10) in || there is a misprint, instead of p GP ^ a one should have /5 GP ^ a on the left 
side.) To see (hi) let us define 

p 7 (x) = 7 2 /(*+ 2 )pTF ( 7 i/(*+2) x ) . (2-22) 



We can write 



with 



e ( 7 ,x) = 7 - s /( s+2 )T/( 7 1 /(«+2) x ) _ jy(x). (2.24) 
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By assumption, |e(7,x)| < o~(7)(l + W(x)) for some 5(7) with lim 7 _ >00 5(7) = 0. Because J" p 7 = 1 
for all 7, we see from Eq. (|2.23[ ) that /u TF (y)7~ s /( s + 2 ) converges to some c as 7 — > 00. Moreover, 
we can conclude that the support of p 7 is for large 7 contained in some bounded set B independent 
of 7. Therefore 



1 = lixn^ j Pi = J ( 87r ) _1 t c - w (*) 



by dominated convergence, so c is equal to the p, of Eq. (1.11). Now 

1 



with 



p 7 (x) = -[A TF -W(x)-6(7,x)]. 



e( 7 ,x) = e( 7 ,x) + /i TF - 7-^+ V F (t)- 



(2.25) 



(2.26) 



(2.27) 



Again |e(7,x)| < <5(7)(1 + W(x)) for some J( 7 ) with lim^,*, $(7) = 0. By Eqs. fll.llD and ( ggg) 
we thus have 



-TFi 



with C = (87T)- 1 sup xeB (l + W(x)) < 00. 

The mean density for the TF theory is defined by 

p 7 = N J ^ F (x) 2 d 2 x. 

For 7 = iV, i.e., g = 1 this is the same as (1.6). It satisfies 
Lemma 2.4 (Bounds on /3 7 ). For some constant C > 



(2.28) 
□ 

(2.29) 



87T7 



7 



(2.30) 



Proof. The upper bound is trivial. Because p 7 , defined in (|2.22j ), converges uniformly to />f F and 
^ TF (^)^/- S /( S + 2 ) _> ^ TF as -y _ > 00, we have the lower bound 



7P 7 



Nfj, (7) 



> 8.7 S/(S+2 V TF (7)-M / (P?I) 2 - 2||p?f 



> C 



for some C > 0. 



(2.31) 
□ 



Remark 2.1. With V asymptotically homogeneous of order s, /i TF (7)7~ s /( s+2 ) converges as 7 — > 
00, i.e. // TF (7) ~ y s /( s + 2 ) for large 7. So the mean TF density for coupling constant g = 1, defined 
in ( |1.6| ), has the asymptotic behavior p ~ jV s /( s+2 ). 
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3 Upper bound to the QM energy 

As in the three dimensional case, cf. Eqs. (3.29) and (3.27) in 0], one has the upper bound 

E^(N,a) J\V4>r\ 2 + V(^ P ) 2 , NJJi^l + lN^M^l f, n 

where we have implicitly used that — A</> GP + V<p^ p > 0, which is justified by Lemma The 
coefficients I, J and K are given by Eqs. (2.4)— (2.10) in [Q]. They depend on the scattering length 
and a parameter b. We choose 7 = iV/| ln(a 2 p)| and b = p~ x l 2 . (Recall that p is short for Pn-) 
With this choice we have (as long as a 2 p < 1) 

4"7T 

J= |ln( a 2p)|' (3 ' 2) 



and the error terms 

<;i> 

NU^Wi.IK (const 



G1 '" 2 ' - < (l + OdM^P)! -1 )) (3.3) 



'7 MOO 

and 

K 2 iV 2 ||^ p ||L < (const.)S GP (l,7)^^ (l + OdMa^l -1 )) , (3-4) 

where we have used Lemma |2.1| . So we have the upper bound 

ffQ M (iV, q) 
£ GP (W,l/|ln(a 2 p)|) 

Now if 7 is fixed as N — > 00 

M GP (7) 



< 1 + O (/^ GP (7)/P) + 0((| In(a^)!- 1 )) . (3.5) 



ln(a 2 p)| iV' 



(3.6) 



If 7 — > 00 with iV we have instead, assuming that the external potential is asymptotically homoge- 
neous of order s, 



M GP ( 7 ) /i TF ( 7 ) / 7 W(s+2) 



/x TF (iV) ViV 



(3.7) 



so m any case 



£ GP (iV,l/|m(a 2 p)|) 
holds as iV — > 00 and a 2 p — > 0. 



<l + |k(a>-)r s ' lw ' (3.8) 



4 Lower bound to the QM energy 

Compared to the treatment of the 3D problem in |^] the new issue here is the TF case, i.e., 
7 = N/\ ln(a 2 /?)| — ► 00, and we discuss this case first. The GP limit with 7 fixed can be treated in 



complete analogy with the 3D case, cf. Remark 4.1 below 
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We introduce again the rescaled jo 7 as in and also 



Note that the scattering length of v is a = a 7 1 '( s + 2 ), Using V > ^ TF ( 7 ) — 87r 7 /9 7 F and (|2.7|) we 
see that 

£Q M (7V,a) > S TF (iV,7/iV)+4^ 7 s /( s + 2 )yp2 +7 -2/( S +2)Q 

-8vriV 7 ^+ 2 )||p 7 -^F|| ooi (4 . 2) 



with 



Q = f| mf 2 / (iV^p + ^^-x.O^^-SvTTpT^xOlvI/l 2 ) . (4.3) 

J 1 1 1 i J \ 3<i J 

Dividing space into boxes a of side length L with Neumann boundary conditions we get 

Q>Y, E hom (n a , L) - 8vr 7 p aimax n a , (4.4) 

a 

where /) ajmax denotes the maximal value of pf\ in the box a, and E hom (n, L) is the energy of a 
homogeneous gas of n bosons in a box of side length L and Neumann boundary conditions. We can 
forget about the boxes where /9 a , m ax = 0, because the energy of particles in these boxes is positive. 
We now want to use the lower bound on E hom given in 0], namely 

E ->,L) > *4]Utkm\ (l-C-IMSW^)r^) • (4.5) 

This bound holds for n > (const.) | ln^n/L 2 )! 1 / 5 and small enough a 2 n/L 2 . Now if the minimum 
in ( |4.4| ) is taken in some box a for some value n a , we have 

E hom (n a + 1, L) - E ho ™(n a , L) > 87r 7Pa , max . (4.6) 

By a computation analogous to the upper bound (see Q) one shows that 

E hom (n + l,L) -E hom (n,L) 

^^ |MgU)l (1 + 0( "" (a ' 2 " /L2rl)) - (4 ' 7> 



Using Lemma 2.4 and the asymptotics of fj, (Remark |2.1[) we see that 



n 2 n n 2 N / N \ 2 /( s + 2 ) „2 n f N\ 2 /( s + 2 ) 



(4.8) 



for some constant C, so ( |4.7| ) reads 

£J hom (n + l,L) -£ hom (n,L) 

n 1 ^ ■ „( i + \Hh/N) 2 ^Lyc)\ \\ 
- 8 ^\H^W\ [ + ° { \H^)\ J J • (49) 

So if L is fixed, our minimizing n a is at least ~ p a ,m&xL 2 N . If JV is large enough and a 2 p is small 
enough, we can thus use (4.5) in (|4.4p to get 



Q > V 4vr [ gf 1 i _ C x _ I . (4.10) 

"-^ U 2 | ln ^^|^ | ln ^^|i/ 5 y |ln(a 2 p)|l 
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Lemma 4.1. For < x, b < 1 we /ioue 

- 2 . b 



x 



lnxl 



| ln6| 



x > 



b 2 



In b\ 



1 + 



1 



(2|ln6|) 2 



Proof. Since lux > — ^x d for all d > we have 



Sr^-2r > ^edx 2 ^-^ >c(d)(6^|ln6|)-VM 
b z \mx\ b b z b 



with 



1 



(2 + (i)(2+d)/(l+d) (2 + d)l/(l+d) 

Choosing d = 1/| In 6| gives the desired result. 
Note that the Lemma above implies for k > 1 



> -1 - -if 
4 



7X&; > 



1 + 



1 



;2j inb|) 5 



I In x| | In fo| | In b 

Applying this with x = a 2 n a /L 2 and b = -/Va 2 /9 Qjmax we get the bound 

Q >-4vriV7^^ max i 2 



1 + 



4 1 ln(a 2 iVp 



\n(a 2 Np 



C 



ln(a 2 p)| 



lnf^) 1 1/5 



for ( 4.1C| ), To estimate the error terms, note that as in ( |4.8| 

■ N \ 2/0+2) 



(4.11) 



(4.12) 

(4.13) 
□ 

(4.14) 



(4.15) 



(4.16) 



so |ln(a 2 iV)| = |ln(a 2 p)| + 0(ln | ln(a 2 p)|) for small a 2 p. Using ||p 7 - PiJ||oo — > (Theorem |2J 
(iii)) and J* p 2 — ► J* (/o^T) 2 as 7 — > oo (which follows from the uniform convergence and boundedness 
of the supports) we get 



hm mi , 
TV^oo E TF (N 



^^ al - (c „ ns ,)(E^-/(«) 2 ) 



(4.17) 



Since this holds for all choices of the boxes a with arbitrary small side length L, and by the 
assumptions on V pfj is continuous and has compact support, we can conclude 



lim inf 



£Q M (iV, a) 



AT^oo ^TF^l/llnCo 2 ^)!) 

in the limit ./V — > 00, a 2 p — ► and N/ \ ln(a 2 /?)[ — > 00. 



> 1 



(4.18) 
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Remark 4.1 (The GP case). In the derivation of the lower bound we have assumed that 7 — » 00 
with N, i.e. N S> [ln(a 2 p)|, which seems natural because otherwise the scattering length would 
have to decrease exponentially with N. However, for fixed 7 one can use the methods of [0] (with 



slight modifications: One uses the 2D bounds on the homogeneous gas and Lemma 4J.) to compute 
a lower bound in terms of the GP energy. The result is 

E QM (N a) 

liminf -^7- /„ ... > 1 (4.19) 
tv^oo E GP (N, 1/1 ln(a 2 /))[) _ v ; 

in the limit N — > 00, a 2 p — > with 7 = N/ \ ln(a 2 p)\ fixed. 



5 The limit theorems 

We have now all the estimates needed for Theorems |1.1| -|1.4|. The upper bound ( |3.8| ) and the lower 
bound (pLl9D prove Theorem |0]. The energy limit Theorem |1.3| for the TF case follows from 



Theorem |2Jj (i) and ( gig) . 

The convergence of the energies implies the convergence of the densities in the usual way by 
variation of the external potential. Replacing V(x) by V(x) + (57 s /( s + 2 )y(^ _1 /( s + 2 )x) for some 



positive Y G and redoing the upper and lower bounds we see that Theorem L3j and Theorem 



2.1 (i) hold with W replaced by W + 5Y. Differentiating with respect to 5 at 5 = yields 



*n ^^P^(7 1/(s+2) x)=^(x) (5.1) 

in the sense of distributions. Since the functions all have norm 1, we can conclude that there is 
even weak L 1 -convergence. 

Remark 5.1 (The 3D case). In the analogues of Theorems [O] and |1.2] were shown for the 
three-dimensional Bose gas. Using the methods developed here one can extend these results to 



analogues of Theorems |l.3| and 1A. In 3D the coupling constant is g = a, so 7 = Na. Moreover, 
the relevant mean 3D density is p 7 ~ ./V(iVa)~ 3 /( s+3 ). 



A Appendix: Scattering length in two dimensions 

Due to the logarithmic behavior of the Green function of the two dimensional Laplacian the def- 
inition of the scattering length is slightly more delicate in two dimensions than in three. For a 
nonnegative potential v(x), depending only on |x| and with finite range Rq, it is naturally defined 
by the following variational principle: 

Theorem A.l. Let R > Rq and consider the functional 

£r[<P] = I {lV0(x)| 2 + ^(x)|^(x)| 2 | d 2 x. (A.l) 

Then, in the subclass of functions such that J(|(/>| 2 + |V0| 2 ) < 00 and </>(x) = 1 for |x| = R, 
there is a unique function 4>q that minimizes £r[4>]. This function is nonnegative and rotationally 
symmetric, and satisfies the equation 

-A0o(x) + ^(x)0 o (x) = O (A.2) 

for |x| < R in the sense of distributions, with boundary condition (po(x) = 1 for |x| = R. 
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For Ro < |x| < R 

o (x) = ln(|x|/a)/ln(i?/a) (A.3) 
for a unique number a called the scattering length. 

For the proof see 0], where generalizations to other dimensions and potentials with a negative 
part are also discussed. Note that the factor ^ in ( A.l ) and (A2) is due to the reduced mass of 
the two body problem. 

If v has infinite range it is easy to extend the definition of the scattering length for nonnegative 
v under the assumption that JZ\ >Rl v(x)d 2 x < oo for some R\. In fact, one may then simply cut 
off the potential at some point Rq > R\ (i.e., set v(x) = for |x| > Rq) and consider the limit of 
the scattering lengths of the cut off potentials as Ro — > oo. See [|| for details. 
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